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Abstract 

A pseudopotential for positronium-atom interaction, based on electron-atom and positron-atom 
phase shifts, is constructed, and the phase shifts for Ps-Kr and Ps-Ar scattering are calculated. 
This approach allows us to extend the Ps-atom cross sections, obtained previously in the impulse 
approximation [Phys. Rev. Lett. 112 , 243201 (2014)], to energies below the Ps ionization thresh¬ 
old. Although experimental data are not available in this low-energy region, our results describe 
well the tendency of the measured cross sections to drop with decreasing velocity at u < 1 a.u. 
Our results show that the effect of the Ps-atom van der Waals interaction is weak compared to 
the polarization interaction in electron-atom and positron-atom scattering. As a result, the Ps 
scattering length for both Ar and Kr is positive, and the Ramsauer-Townsend minimum is not 
observed for Ps scattering from these targets. This makes Ps scattering quite different from elec¬ 
tron scattering in the low-energy region, in contrast to the intermediate energy range from the Ps 
ionization threshold up to u ~ 2 a.u., where the two are similar. 

PACS numbers: 34.80.-i, 36.10.Dr 
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I. INTRODUCTION 


Recently observed similarities between the positronium scattering and the electron scat¬ 
tering from a number of atoms and molecules m in the intermediate energy range suggest 
that both processes are largely controlled by the same interactions. When plotted as a 
function of the projectile velocity, the electron and Ps cross sections are very close and 
even show similar resonancelike features. This similarity was explained recently |1] using 
the impulse approximation. It was shown that above the Ps ionization threshold, the Ps-A 
interaction is mainly controlled by the e~-A scattering amplitude, and the e~-A exchange 
contributes mostly to this amplitude in the intermediate energy range. It is clear, however, 
that at higher energies, the role of the exchange interaction becomes less significant, and 
the similarity between Ps-A and e~-A scattering should gradually disappear. On the other 
hand, at lower energies, long-range interactions between the projectile and the target play a 
signihcant role, and they are substantially different for e~-A and Ps-y4 interactions. In the 
former case it is the polarization potential decreasing as r~^ at large distances r, and in the 
latter case the van der Waals interaction decreasing as r“®. 

The impulse approximation for Ps-A scattering |1] produces very large cross sections 
below the ionization threshold because of the dominance of the e^-A scattering amplitude. 
This growth of the Ps cross section is unphysical, since the large e'^-A amplitude is due to 
the effects of positron-atom polarization and virtual Ps formation, both of which are absent 
in Ps-A scattering. Since the impulse approximation breaks down at energies below the 
ionization threshold, alternative methods and approximations, such as close-coupling and 
static-exchange, should be used in this energy range. Blackwood et al. |5] performed close- 
coupling calculations of Ps scattering from noble-gas atoms. These calculations allowed for 
the distortion and break-up of Ps, but kept the target “frozen”, i.e., they neglected any 
excitations of the target. Virtual target excitations are known to be very important in low- 
energy electron- and positron-atom collisions, where they can be described in terms of the 
polarization interaction. For low-energy Ps-atom collisions, they give rise to the van der 
Waals interaction, which can be incorporated by extending the close-coupling calculations 
to include the virtual excitations of the target. Such calculations have been performed for 
Ps collisions with the hydrogen atom laiii and would be an ultimate goal in the problem of 
Ps-atom collisions. However, they are very challenging computationally for complex atoms. 
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In this paper we develop an alternative low-energy method based on the use of the electron 
and positron scattering phase shifts, similar to the impulse approximation. It involves 
constructing model potentials that reproduce these phase shifts, and then adding them 
to describe the Ps-atom interaction. While this procedure is straightforward for positron 
scattering, the situation with electrons is more complicated. Due to the Pauli exclusion 
principle, the effective potential for the electron depends on its orbital angular momentum 
Z, i.e., it becomes a pseudopotential |8]. When such a pseudopotential is averaged over the 
electron density distribution in the Ps atom, it becomes a nonlocal operator. 

Another difficulty is related to inclusion of the long-range interaction. Accurate low- 
energy electron and positron scattering phase shifts contain contributions of the atomic 
polarization potential. This potential is attractive for both electrons and positrons and 
behaves as —a/2r^ at large distances, where a is the atomic dipole polarizability. (We use 
atomic units throughout.) An effective Ps-A potential including such contributions would 
behave as —a/r^, which is physically incorrect, as the dominant long-range Ps-A interaction 
is the van der Waals interaction —Cq/R^. The latter potential results from the many- 
body Ps-atom, rather than single-particle (i.e., electron-atom or positron-atom) dynamics. 
It can be obtained by including the two-body polarization potential (see, e.g., 0) in the 
Hamiltonian, which gives the total polarization interaction at large distances as 


^poi(re,rp) 


are - Tp 

^ 3^3 
e p 


( 1 ) 


where rg and are the electron and positron position vectors, respectively, relative to the 
target. Averaging of this potential over the electron and positron density distribution in Ps 
does lead to an effective van der Waals interaction HD]. 

Alternatively, one can construct the positron-atom and electron-atom pseudopotentials 
using static (static-exchange) phase shifts for the positron (electron), i.e., neglecting the 
polarization. The Ps-atom van der Waals interaction can then be added, e.g., in the form 


C'e 




( 2 ) 


where Cq is the van der Waals constant and Rc is a cutoff radius. The Cq values for Ps- 
atom pairs are known, e.g., from the London formula which gives Cq = 152 a.u. for 
Kr and Cq = 104.5 a.u. for Ar. (These values are close to the estimates obtained in Ref. 
[12].) In contrast, the cutoff parameter Re cannot be determined rigorously, but the phase 
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shifts and cross sections are sensitive to its choice. A similar problem is enconntered when 
nsing the polarization interaction Q, which also reqnires a cnt-off at small distances. In 
the present calculations the radius Rc is set by requiring that the cross sections given by 
the pseudopotential method merge smoothly with the elastic cross section calculated in the 
impulse approximation above the Ps ionization threshold |1]. 

The rest of the article is organized as follows. First we discuss the construction of the 
pseudopotentials for e+ and e~ scattering from the static (static-exchange) phase shifts, and 
the derivation of the pseudopotential for Ps-atom scattering. We then present the results 
for Ps-Kr and Ps-Ar scattering and discuss the low-energy behavior of the cross sections. 
For both atoms, the scattering length is positive, which implies effective repulsion at low 
energies and rules out the existence of the Ramsauer-Townsend minimum. 


II. THEORY 


A. Pseudopotentials 


We choose the positron-atom pseudopotential in the form 


V,{t) = 


P „—a„r 


( 3 ) 


which represents the static e~^-A repulsion, and where Zp and ap are htting parameters. They 
are obtained by htting the s-, p- and d-wave scattering phase shifts in the potential ([^ to the 
positron scattering phase shifts in the static held of the ground-state atom calculated in the 


Hartee-Fock approximation (see Sec. IIB). The parameter Zp plays the role of an ehective 
nuclear charge. It can be diherent from the actual nuclear charge, since a low-energy positron 
does not penetrate deep into the atom. 

The ehective static-exchange potential for the electron is chosen as 




( 4 ) 


where the second term represents repulsion due to the Pauli exclusion principle. This ehect 
depends on the orbitals occupied in the atomic ground state, hence it is /-dependent. We 
also found that in general the ab initio static-exchange (Hartree-Fock) phase shifts cannot 
be reproduced using Ze = Zp. Therefore, we regard all parameters in Eq. (|^ as /-dependent. 
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Formally this means that the effective electron-atom potential is a nonlocal operator with 
the kernel 

K(r, r') = \/{r - r') (5) 

Im 

where Vi{r) are potentials given by Eq. Q. It is convenient to rewrite this expression as 
l/e(r,r') = -Vp{r)6{r - r') + ^6{r - r')J2Mr)yi*m{r)Yim{r'), (6) 

Im 

where 

Vi{r) = Vi{r) + Vp{r). (7) 

Since the “direct” part of the potential Vi{r) is close in magnitude, but opposite in sign to 
Vp{r), vi{r) represents mainly the exchange interaction between the electron and the atom. 
The Ps-atom pseudopotential can be now written as 

Im 

( 8 ) 

In the static approximation we average this potential over the electron and positron density 
distribution in Ps given by |<F(p)p, where $(p) is the Ps ground-state wave function and p 
is the relative coordinate. The relations between re, fp and p are 


l/ps(re, r;, rp, r') = Vp{rp)S{rp - - Vp{re)S{re - r'J + — 


r^ = R + p/2, rp = R-p/2, 


where R is the position of the Ps center of mass relative to the target. 

The average of the local part of the pseudopotential (|^ reduces to the integral 

j [Ep(R-p/2)-Ep(R + p/2)]|$(p)|>, (9) 

which vanishes because the integrand is parity odd. This corresponds to a well-known fact 
that the static potential for the Ps-A interaction is zero. The remaining nonlocal part in 
Eq. ([^ contains a strong repulsive core, and to make the calculations more tractable, it is 
convenient to represent vi{re) as 


Vl{re) = ViociVe) +Mz(re), 

where the /-independent part fioc(’"e) contains the major repulsive contribution, and ui{re) 
accounts for the remaining /-dependent part. The averaging procedure is then reduced to 


5 


averaging of the Ps pseudopotential 


Vps(re, r;,) = !ii„c(r,)(5(r, - r'J + -^i(re - r'J 


( 10 ) 


Im 


The average of the local part of this potential gives the local Ps-atom potential 

Wav(^) = [ nioc(R. +p/2)|<h(p)|^cip, 


( 11 ) 


while averaging the nonlocal part [second term in Eq. (10)] gives a nonlocal contribution to 
the Ps-A interaction, 


V'(R,R') = ^ f ii(r-r>,(r)F,;„(f)yi„(f')|4>(p)|= 

Im 


dp 


( 12 ) 


where r = R + p/2 and r' = R' + p/2. 

In performing these integrations, it is more convenient to switch to the integration variable 
r. Since the result depends only on the absolute magnitude of the vectors R and R', and 


on the angle between them, the integration in Eq. (12) can be performed in the coordinate 
system with the polar axis along the vector 

s = R' - R. 


We then have 


V (R, RQ = 8 ^ [ ^6{r - r')ui{r)Pi{cos6rr')\^{2{r - K))\'^dr, 


dvr / P 


(13) 


where r' = r + s and 6rr' is the angle between r and r'. 

Integration over the polar angle 6 in this coordinate system eliminates the ^-function 
(which ensures r = r') and Exes the angles: 


s „ r + s cos 9 ^ s" 


cos 9 = -, cos 9rr/ = 

2r |r + s 


= 1 - 


2r2' 


It also introduces a factor 1/s since 


d{r — |r + s|) 
dcos9 

for cos 9 = —sj2r. 

The Ps ground-state density is expanded in spherical harmonics as 

oo 

|$(2|r-R|)|2 = —F,(r, R)(2/'+ l)P,(cos0rR), 


(14) 


l'=0 
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where the expression for Fi (r, R) is given in Appendix (see also Appendix B in Ref. m)- 
The Legendre polynomial Pii {cos 6 ^ 11 ) is the only part of expansion (14) which depends 


on the azimuthal angle 0 in the integrand of Eq. (13) (^rr' does not depend on 0). Therefore 
we can perform integration over 0 as 

r 2 TT 


Pi'{cos 6 r-R)d(j) = 27 tP i'{cos 6 )Pi> {cos 6 r), 


where 


cos6r = 


P! cos Q — R 


and 0 is the angle between R and R'. 

Thus, we obtain the nonlocal part of the Ps-atom potential as 

P(R,R') = ^J2i‘^l + l){2l' + l)Pi'{coseR) r Pi (l - Pi' (-^) Fi'{r,R)ui{r)dr. 
IV ^0 V / 

( 16 ) 

It is convenient to expand this expression in Legendre polynomials, 

V(R,R') = — —,^l4(R,R')-Pi(cose), 


PR' ^ dvr 

L=0 


where 


n'K 

Vl{R,P') = 2ttRR' / l/(R,R')PL(cos0)sin0d0. 

Jo 

Similarly, for the local part of the interaction potential, given by Eq. 01 , we obtain 


nav(R)= 4 / vioc{r)Fo{r,R)r‘^dr. 

Jo 

Substitution of the local and nonlocal potentials in the Schrodinger equation for the Ps-A 
system yields a set of radial equations: 


1 d^h 

2 m dR^ 


+ 


E -Vav{R) 


L{L + 1) 

2 mR^ 


fdP) - / Vl{R, P!)fL{P!)dP! = 0, (16) 


where m = 2 a.u. is the Ps mass, and fL{R) is the radial part of the Ps center-of-mass wave 
function for the orbital angular momentum L. 


The sums in Eq. (15) converge fast. With a proper choice of nioc(’"e), the sum over I can 
be truncated at /max = 2. Increasing /max to 4 has almost no effect on the phase shifts in the 
velocity range up to 2 a.u. The sum over V converges if > 6. 
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B. Fitting parameters 


The values of the parameters of the pseudopotentials in Eqs. 0 and 0 , were determined 
by fitting the positron-atom and electron-atom scattering phase shifts obtained in the static 
potential of the atom calculated in the Hartree-Fock approximation m- 


As mentioned in Sec. IIA for the positron, the same pseudopotential can be used for all 
partial waves. The static-held scattering phase shifts for the positron on Ar are shown in 
Fig.|Il and the corresponding sets of parameters for Ar and Kr are given in Table For 
positron velocities up to 2 a.u., the pseudopotential phase shifts are within 1% of the actual 
static-held phase shifts and are indistinguishable from them on the scale of the plot. 



FIG. 1: Positron scattering phase shifts for Ar obtained using the static potential of the Hartree- 
Fock atom. 

For the electron, the pseudopotential is optimized separately for each partial wave, by 
htting the Hartree-Fock scattering phase shifts. The corresponding sets of parameters for Ar 
and Kr are given in Table [T| and the phase shifts for Ar are shown in Fig. for I = 0-2. The 
pseudopotential for the s wave gives the phase shifts that are almost indistinguishable from 





TABLE I: Parameters of the positron and electron pseudopotentials, given 

Ar and Kr. 

by Eqs. ([^ and (Q, for 

System 

1 

Z 


B 

n 

(3 

e''‘-Ar 

0-4 

18.06 

1.95 

- 

- 

- 

e“-Ar 

0 

4.297 

0.618 

14.72 

3 

0.248 


1 

10.0 

2.368 

15.24 

3 

0.504 


2 

9.780 

1.230 

0 

- 

- 


3 

12.48 

1.544 

0 

- 

- 


4 

15.13 

1.714 

0 

- 

- 

e^-Kr 

0-4 

20.79 

1.760 

- 

- 

- 

e“-Kr 

0 

20.79 

1.760 

56.84 

6 

0 


1 

20.79 

1.760 

97.80 

6 

0 


2 

18.25 

1.317 

0 

- 

- 


3 

14.80 

1.409 

0 

- 

- 


4 

16.98 

1.544 

0 

- 

- 


the Hartree-Fock, and the htted phase shifts for the p and d waves are also qnite accnrate. 
The behavior of the positron and electron phase shifts and the qnality of the psendopotential 
hts for Kr are similar to those shown for Ar. For both atoms, the d-wave displays a broad 
resonance at the electron velocity n ~ 1 a.n. (see Fig. |^. 


III. RESULTS AND DISCUSSION 


In this section, we present the Ps-atom scattering phase 
which contains the static Ps psendopotential, and with the 
interaction (|^. We also nse these phase shifts to compnte 
cross section, and compare onr resnlts with experiment and 


shifts calcnlated from Eq. (16) 
inclnsion of the van der Waals 
the Ps-atom elastic scattering 
other theories. 
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electron velocity (a.u.) 


FIG. 2: Electron scattering phase shifts (modulo vr) for Ar obtained using the static-exchange (i.e., 
Hartee-Fock) atomic potential (solid black lines), and the pseudopotential Q (dashed red lines), 
with parameters listed in Table |I] 

A. Ps-Kr scattering phase shifts 


Figure shows the phase shifts for Ps with the orbital angular momentum L = 0, 1, 
and 2, scattered from Kr. Adding the van der Waals interaction (with Rc = 3.0 a.u., see 


Sec. Ill C) produces relatively small positive corrections for the s and p waves, while for the 
d wave the correction is relatively large. In particular, due to the van der Waals potential, 


the phase shift acquires a characteristic rise, 82 oc at low Ps momenta k [see Eq. (B6) . 
Our calculations show that for L = 0 and 1, the nonlocal effects related to the second 


term on the right-hand side of Eq. (10) are small, but for L = 2 they are significant. This 


can be seen from Fig. which compares the full (nonlocal) d-wave phase shift with those 
obtained using two choices of the local potential, namely, uioc(’") = vi{r) and uioc(’") = V 2 {r), 
where vi{r) is defined by Eq. ([^. The choice of uioc = V 2 leads to a shape resonance due to 
the combination of the van der Waals interaction and the centrifugal barrier. The resonance 
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FIG. 3: (Color online) Ps-Kr scattering phase shifts for L = 0, 1 and 2. For each partial wave, 
the lower curve (solid black) is for the static pseudopotential, and the upper curve (dashed red) 
includes the effect of the van der Waals interaction with Cq = 152 a.u. and Rc = 3.0 a.u. 


becomes very pronounced at smaller Rc. However, for uioc = vi, the resonance is not visible, 
and the result obtained with the full potential (i.e., including the local and nonlocal terms) 
conhrms that the nonlocal effects suppress the resonance behavior of the Ps d wave. This 
suppression leads to a better agreement with the static-exchange calculations of Blackwood 
et al. [S] (see below). A trace of this resonance is a weak minimum in the L = 2 partial 


cross section (see Sec. IIIC below). 

Analysis of the Ps s-wave phase shift at low momenta (^o — —Ak) yields the value of the 
scattering length A. In the static approximation we hnd A = 3.32 a.u., and when the van 
der Waals interaction is included we obtain A = 2.35 a.u. (for Rc = 3.0), or A = 2.50 a.u. 
(for Rc = 3.5). These values can be compared with those of Mitroy and Bromley [13], i.e., 
A = 3.18 a.u. in the static approximation, and A = 1.98 a.u. with the van der Waals 
interaction included. Note that because of the uncertainty in the van der Waals interaction 
effect in Ref. na, the corresponding scattering length varies between 1.22 and 2.26 a.u., 
with A = 1.98 a.u. being their best prediction. 
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FIG. 4: (Color online) Ps-Kr scattering phase shifts for L = 2. Shown are the results obtained with 
the full pseudopotential (i.e., with nonlocal effects) and its local approximations uioc(?') = vi{r) for 
I = 1 and 2. The van der Waals interaction is included in all cases. Solid black lines: Rc = 3.0 a.u., 
red dashed line: Rc = 3.5 a.u. 

B. Effect of the van der Waals interaction on the scattering length 


As seen from the s-wave phase shift in Fig. the effect of the van der Waals interaction is 
not very significant. It does not change the sign of the scattering length (the way atomic po¬ 
larization does for the electron and positron scattering), and there is no Ramsauer-Townsend 
minimum in the cross section. To understand this qualitatively and semiquantitatively, con¬ 
sider a model potential with a hard repulsive core of radius Rq and a van der Waals “tail”, 




(17) 


-f-oo, R <Z Rq 

— CqIR^, R > Rq 

Here the repulsive core mimics the exchange interaction between the electron and the atom. 
The s-wave radial Schrodinger equation for A: = 0 at R > i?o is 

dVo 27716-6 

- 0 - 


Its solution (up to a normalization constant) is 


/o(R) = [J-i/4{xo)Ji/a{x) - Ji/a{xo)J-i/4{x)] , (18) 
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where 


X = 






Xo = 


\fmC^ 


Rl 


(19) 


and Jy is the Bessel function. The scattering length, obtained from the asymptotic behavior 


of Eq. (18) at i? —)■ oo, /o(i?) (x R — A, is 
mCe r(3/4) J_i/4(a;o) 


A = 


r(5/4) Ji/4(a:o) 


2 a; 


= 1 - - 


22 x^ 


844 Xg 


15 1575 482625 


+ ... 


( 20 ) 


Here the factor (mC'6/8)^/‘^r(3/4)/r(5/4) is similar to the mean atom-atom scattering length 
determined by the long-range part of the van der Waals interaction (see Ref. dS]). 

Figure]^ shows the scattering length (20) as a function of Rq for three values of the van 
der Waals constant, Ce = 104.5, 152, and 234 a.u., corresponding to Ps-Ar, Ps-Kr, and 
Ps-Xe interactions. The scattering length becomes negative only for unrealistically small 
values of Rq. As an estimate, we can assume that Rq equals the scattering length in the 


static approximation. For Kr, this gives Rq = 3.32 a.u. Using this value in Eq. (20), we see 
that the van der Waals force reduces the scattering length to A = 2.67 a.u., in good accord 


with the calculations (Sec. Ill A). Therefore, the effect of the van der Waals interaction is 
not as drastic as the effect of the polarization interaction in e“-Kr (or e^-Kr) scattering, 
where it makes the scattering length negative. Obviously this is due to the shorter range of 
the van der Waals force as compared to the polarization force. 

Alternatively, one can estimate Rq using the mean radius of the outer atomic orbital (r), 
e.g., as 

Ro = ^{r), ( 21 ) 

where 7 ~ 1 is a dimensionless factor. For Kr, (r) = 1.95 a.u. nzi, and in order to obtain 
the correct scattering length, A = 2.35 a.u. (for Rc = 3.0), we should choose 7 = 1.61, 
which leads to Rq = 3.14 a.u., close to our previous estimate of Rq = 3.32 a.u. 

We can use this simple model to estimate the Ps scattering lengths for other atoms. 
Table [H] lists the values of (r) from Ref. m and the corresponding values of the model 


scattering length (20) for Ar, Kr and Xe, obtained using 7 = 1.61. Also shown are the 
results of the scattering calculations of this work and of Mitroy et al. urn- 

Although the model result for Xe looks somewhat overestimated, the model correctly 
predicts the trend of the scattering length A to increase with the atomic number Z. Although 
the van der Waals interaction (which makes A smaller) grows with Z, the increase of the 
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FIG. 5: Dependence of the scattering length, given by Eq. (20), for the van der Waals potential 
© on the radius of the repulsive core. The three curves correspond to Cq = 104.5, 152, and 234 
a.u., for Ps-Ar, Ps-Kr, and Ps-Xe, respectively. Open circles show the estimates of A from Table 

m 


TABLE II: Mean atomic radii (r), core radii Rq, and Ps-atom scattering lengths A for Ar, Kr and 
Xe. All values are in a.u. 


System 


(r) 

Rq 



A“ 

Ar 

104.5 

1.66 

2.67 

1.73 

2.14-2.33 

1.30-1.98 

Kr 

152 

1.95 

3.14 

2.35 

2.35-2.50 

1.22-2.26 

Xe 

234 

2.39 

3.85 

3.23 

- 

1.50-2.60 


“Scattering length from Eq. (20), obtained using Eq. (21) with 7 = 1.61. 
^Present scattering calculations. 


“Values obtained by Mitroy et al. [siiin]. 


atomic radius, i.e., the effect of the static repulsion, is stronger. [This can be seen from the 
expression for xq, given by Eq. (19), which is proportional to but inversely proportional 
to Rq] This is opposite to what is observed in electron-atom scattering where the effect of 
the polarization attraction takes over the effect of the core radius, and the scattering length 
decreases with the growth of Z. (The analog of the parameter xq in this case is proportional 
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to -y/a/Ro, see Ref. [T6].) 

These observations have an important consequence for the comparison of e~-A scattering 
with Ps-y4 scattering. The observed similarities [T] at energies above the Ps ionization 
threshold are explained in terms of the impulse approximation |1]. On the other hand, in the 
region below the ionization threshold, where the impulse approximation fails, no similarity 
exists. In this energy range scattering is controlled by different long-range interactions, 
the strong polarization interaction for the electrons, and the relatively weak van der Waals 
interaction in the case of Ps. 

C. Ps-Kr scattering cross section 

Figure shows the elastic Ps-Kr cross sections in the velocity range from threshold to 
2 a.u. Higher partial waves (L > 4), up to L = 10, were included by solving the radial 
equation in the local approximation with inclusion of the van der Waals interaction. The 
van der Waals interaction was included with two choices of the cut-off parameter: Rc = 3.0 
and 3.5 a.u. The figure also shows the cross section obtained by Blackwood et al. |5] using 
the static-exchange approximation, i.e., without inclusion of virtual excitations in the target 
or projectile. These results can be compared with our static calculations. Figure shows 
that the two theories are very close at low velocities, but the present static-held cross section 
decreases faster with the increasing velocity, compared to that of Blackwood et al. However, 
after inclusion of the van der Waals interaction, our cross section increases signihcantly at 
V > 0.6 a.u., and merges with the result of the impulse approximation for Rc = 3.0 a.u. We 
therefore choose this value of Rc for comparison with the experiment (see below). 

At low velocities (below the Ps ionization threshold, v = 0.5 a.u.), the van der Waals 
interaction leads to a signihcant reduction of the elastic cross section. This effect is due to 
the decrease in the scattering length and a general increase of the low-/ phase shifts (i.e., 
decrease in absolute magnitude), as seen in Fig. The van der Waals interaction also 
results in two features in the cross section: a local maximum at n = 0.22 a.u. and a local 
minimum at n = 0.67 a.u. (for Rc = 3.0 a.u.). The former is due to the s-wave contribution 
and is caused by the long-range attractive interaction. To understand this, it is useful to 
discuss a similar effect in the low-energy electron-atom scattering which is controlled by the 
polarization interaction. According to the modihed effective-range theory of O’Malley et al. 
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FIG. 6: (Color online) Ps-Kr elastic scattering cross sections. Solid black curve “lA” is the result of 
the impulse approximation [1]; solid black curve “static” is the present static-field calculation (i.e., 
Cg = 0); dashed red curve is the present calculation with the static and van der Waals interaction 
{Cq = 152 a.u., Rc = 3.0 a.u.); dotted blue curve is the same for Rc = 3.5 a.u.; dot-dashed magenta 
curve is the static-exchange calculations of Blackwood et al. [3]. 


Is] , at low energies the s-wave phase shift behaves as 


tan 5q = —Ak — vr/c^a/S -|- 0 {k^ In k), 


( 22 ) 


where a is the atomic polarizability. Here the characteristic quadratic part of the polarization 


contribution to the phase shift [second term in Eq. (22)] is negative (although the total 
contribution of the attractive polarization potential is positive). As a result, for a negative 
scattering length A, the phase shift passes through zero at small A; > 0 , leading to the 
Ramsauer-Townsend effect. In contrast, for A > 0 the phase shift decreases faster than 
linear, which gives rise to a maximum in the partial cross section. For example, a maximum 
is observed in electron scattering from Ne [12], for which the scattering length is small and 
positive. (Note that although the maximum in the total cross section for Ne is observed at 
about E = 25 eV, the s-wave contribution peaks at E = 6.7 eV, still quite a large energy 
compared to the position of a typical Ramsauer-Townsend minimum.) One could call this 
phenomenon the “anti-Ramsauer” effect, although we are not aware of the use of such term 
in the literature. 
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A similar situation occurs in Ps-atom scattering, although now the additional contribution 
to the phase shift comes from the van der Waals interaction. As shown in Appendix the 
modified effective-range expansion of the s-wave scattering phase shift reads as 


tan(5o = —Ak — Bk^ -\- 2mCQ7ik'^/15 + 0{k^), 


(23) 


where B = and rg is the effective range [SO]- The coefficient B depends on both the 

short-range and the long-range (van der Waals) interactions. If B > 0, the major correction 
to the —Ak behavior is negative, since the k‘^ term is relatively small at low energies. Our 


calculations show that this is indeed the case, although the expansion (23) is valid only at 


very low energies. Hence, a weak anti-Ramsauer effect is observed. Naturally, the effect is 
small compared to that observed in electron-atom scattering because of the relative weakness 
of the van der Waals interaction compared to the polarization interaction. 

The minimum at n = 0.67 a.u. observed in Fig. [^for Rc = 3.0 a.u. is due to the d-wave 
shape resonance, which is quite pronounced in the e“-Kr scattering m However, in Ps-Kr 
scattering, this resonance is suppressed, as discussed in Sec. |IH A| and seen in Figs. and 
1^ Due to a strong background contribution, this resonance appears as a window. 

To compare with the experimental Ps-atom total scattering cross section, the Ps ion¬ 
ization cross section should be added to the elastic cross section. Indeed, Ps ionization 
contributes significantly at velocities v > 0.5 a.u. As in the impulse-approximation calcula¬ 
tions |1], the ionization cross sections are taken from Ref. |22]. Figure shows the elastic 
and total cross section computed in the present work together with the results of the impulse 
approximation and experiment. Although the experimental data [1] are not available at low 
velocities, the data point at n = 0.63 a.u. indicates that the cross section should slightly 
decrease towards lower velocities. This trend is confirmed by our results. The peaking of the 
experimental cross section at n ~ 0.9 a.u. is also in agreement with our results. In addition, 
our calculation predicts a weak local maximum at n = 0.23 a.u. and a local minimum at 
V = 0.56 a.u. Both of these predictions call for experimental verification. 


D. Ps-Ar scattering 

Figure 1^ shows the s-, p- and d-wave phase shifts for Ps-Ar scattering. All of the main 
features here are the same as in the Ps-Kr scattering. However, tan do follows the linear 
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FIG. 7: (Color online) Ps-Kr total scattering cross sections. Solid black curve is the total cross 
section in the impulse approximation [3]; dashed red curves present elastic and total cross sections 
obtained with Cq = 152 a.u., Rc = 3.0 a.u. Solid squares are experimental data from Ref. [T]. 

—Ak behavior much more closely. This results in a decrease of the s-wave cross section, 
according to 

where tq is, in fact, quite small. At the same time the p-wave contribution increases rather 
sharply from threshold, leading to a local maximum in the total cross section at n = 0.22 a.u. 
in the static approximation {Cq = 0), and an even more pronounced local maximum at 
V = 0.26 a.u. when the van der Waals interaction is included. The total elastic cross 
sections are shown in Fig. 

The local maximum in the total cross section is controlled by the local part of the Ps-Ar 
pseudopotential for L = 0 and L = 1. Since this potential is very sensitive to the partial 
cancellation of the attractive part due to the e“-Ar interaction and the repulsive part due 
to the e’''-Ar interaction, the position and the magnitude of the maximum is subject to 
uncertainties. It is possible that effects not included in the present calculations, e.g., short- 
range correlations, can change the position and shape of the maximum or even eliminate 
it completely. We note that the close-coupling calculations of Blackwood et al. [5], which 
allow for the virtual excitations of the Ps, but keep the target atom frozen, do not exhibit 
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FIG. 8: (Color online) Ps-Ar scattering phase shifts. For each L, the lower curve (solid black) was 
calculated using the static-field pseudopotential, while the upper (dashed red) also includes the 
effect of the van der Waals interaction with Cq = 104.52 a.u. and Rc = 3.0 a.u. 


a low-energy maximum at all. However, they do show a local maximum at u = 1.1 a.u., 
which in our calculation is due to the d-wave contribution. Overall, Fig. shows that there 
is a reasonable agreement with calculations [5] for v = 0.4-1 .2 a.u.. One can also see that 
the effect of the van der Waals potential for Ar is smaller than it is for Kr (Fig. |^, and the 
cross section is not as sensitive to the choice of Rc (which can be related to the smaller Cq 
value for Ar). For a more detailed comparison between the present calculation and that of 
Blackwood et al. |S], information on partial cross sections would be needed. 

The Ps-Ar scattering length in the static approximation is A = 3.19 a.u., whereas the 
value calculated by Mitroy and Ivanov [9] is 2.85 a.u. After adding the van der Waals 
interaction we obtain A = 2.33 a.u. for Rc = 3.0 a.u., and A = 2.14 a.u. for Rc = 2.5 a.u. 
The corresponding value obtained by Mitroy and Ivanov varies between 1.30 and 1.98 a.u. 


Comparing with the Ps scattering lengths obtained for Kr (Sec. Ill A), we see that the values 
for both atoms are quite close. While the Ar-Kr system has a smaller Cg value, the Ar atom 
has a smaller radius, and the two effects largely compensate for each other (see Fig. [^and 
Table 0. 

As seen in Fig. the choice of Rc = 3.0 a.u. matches better with the impulse ap- 
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FIG. 9: (Color online) Ps-Ar elastic scattering cross sections. Solid black curve “lA” is the result of 
the impulse approximation [3]; solid black curve “static” is the present static calculation (Cg = 0); 
dashed red curve is the calculation with the static and van der Waals potential (Cg = 104.5 a.u., 
Rc = 2.5 a.u.); dotted blue curve is the same for Rc = 3.0 a.u.; dot-dashed magenta curve is the 
close-coupling calculations of Blackwood et al. [5] . 

proximation results at higher velocity. However, the smaller value of Rc = 2.5 a.u. leads 
to a better agreement with the experiment. Figure shows the total cross sections for 
Ps-Ar scattering obtained for each of these cut-off radii. Whereas both theoretical curves 
describe well the overall behavior of the measured cross sections [23] at u > 0.5 a.u., the 
cross section obtained with Rc = 2.5 a.u. agrees better with the measured absolute values. 
This value also appears to be more physical, as it is smaller than the optimal cut-off radius 
of Rc = 3.0 a.u. found for Kr. 

IV. CONCLUSION 

The pseudopotential model developed in the present paper describes Ps scattering from 
noble-gas atoms (Ar and Kr) at energies below the ionization threshold and matches the 
impulse-approximation results above the ionization threshold. Although experimental data 
are not available in the low-energy region, our results describe well the trend seen in the 
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FIG. 10; (Color online) Ps-Ar total scattering cross sections. Solid black curve “lA” is the result 
of the impulse approximation [1]; dashed red curve is the present elastic cross section obtained 
with Ce = 104.5 a.u., Rc = 2.5 a.u., augmented by the ionization cross sections of Starrett et al. 
[22]; dotted blue curve is the same for Rc = 3.0 a.u. Open circles: experiment [23]; solid squares: 
experiment [1]. 

experimental cross sections to drop with decreasing velocity below n 1 a.u. In addition, 
our calculations predict zero-energy cross sections (or scattering lengths) which are in accord 
with stochastic variational calculations mm- 

Analysis of the scattering phase shifts shows that the static Ps-atom interaction is repul¬ 
sive. This repulsion arises from the electron Pauli exclusion from closed-shell atoms (while 
the pure electrostatic interaction is zero for the truly neutral Ps atom). The phase shifts 
also indicate that the role of correlations represented by the van der Waals interaction at 
low energies, is relatively small. 

Because of the relative weakness of the van der Waals interaction compared to the po¬ 
larization interaction in electron-atom scattering, the scattering lengths for both Ar and 
Kr are positive, and the Ramsauer-Townsend minimum is not observed for these targets. 
The overall picture of Ps-A scattering is quite different from the e~-A scattering in the 
low-energy region. This is in stark contrast to the intermediate energy range from the Ps 
ionization threshold up to u ~ 2 a.u. Here the Ps-A scattering is mostly controlled by the 
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electron-atom exchange, which makes its cross section very similar to that for e~-A scatter¬ 
ing. In the low-energy region, where the long-range interaction is important (especially for 
the electrons), this similarity disappears. 

Although the van der Waals interaction in Ps-A scattering does not produce the 
Ramsauer-Townsend minimum, it can lead to more subtle features in the cross sections, 
such as low-energy maxima. However these features are subject to uncertainties because of 
a delicate balance between the repulsive and attractive components of the Ps-A interaction. 
They can also be affected by other effects such as short-range correlations, which are not 
included explicitly in the present calculation. The uncertainty can be resolved by perform¬ 
ing accurate measurements of Ps-A scattering at low energies, and by new fully correlated 
calculations for this interesting and challenging system. 
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Appendix A: Expansion of the Ps ground state density 

We start with a known expression for the free-particle Green’s function (see, e.g., m), 

„ikri2 


ri2 


ik'^{2l + l)h[^\x>)ji{x<)Pi{cos9i2), 


1=0 


where x> = fcr>, x< = /cr<, and ji, are the spherical Bessel and Hankel functions. 

To switch to the decaying exponent, we make a substitution k = in and arrive at the 
following expansion 

e 


—12 ^ J2{21 + l)FKri,r2)Pz(cos0i2) (Al) 


1=0 


where 


Pi(^iA2) = ^ Khp^(mr>)j)(mr<) 


It is convenient now to introduce the following real functions: 

hi{x) = -i^h\^\ix), ji{x) = i^jiiix), 
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which are related to the modified Bessel function Kij^ii 2 {x) and //+i/ 2 (a;). Explicit expres¬ 
sions for the first few of these functions are 

ho{x) = -, hi{x) = -(l-h-J, h2{x) = -(iH- 

X X \ X J X \ X x^ 

sinh X - If sinh x \ - sinh x f 3 \ 3 cosh x 

jo{x) = -, ji{x) = - -coshx , j 2 {x) = - I 1 + — 


X ' ' X \ X 

The recurrence relations for these functions are 


X 


x^ 


x^ 


Zi+i{x) = ^ ^ Zi{x) -f Zi_i{x), 


X 


dzfx) I + 1 


zfx) - Zi-i{x), 


dx X 

where zi stands for either ji or hi. Using these, we obtain 


Fo{xi,X 2) = -ho{x>)jo{x<) + x>hi{x>)jo{x<) + a;<ho(a:>)ii(a:<), 


and 


Fi{xi,X 2 ) = {-iy[{2l l)hi{xy)ji{x^) +a;>h;_i(a;>)j/(a;<) -1-a;<hi(a;>)j/_i(a;<)], I > 1. 


In practice, expansion (Al) converges well by summing up to /max = 30 if, for high /, 


asymptotic expansions for hi and ji are used. 


Appendix B: Modified effective range expansion for the van der Waals potential 


We are interested in the behavior of the s-wave scattering phase shift (5o. According to the 
effective range theory for short-range potentials, it is given by the effective-range expansion 


fc cot ^0 = “4 + + 0{k^) 

/]. z 


(Bl) 


where A is the scattering length and vq is the effective range. Note that in the presence of 
a weakly bound state, rg > 0 [ 20 ], but generally this is not true. For small momenta and 


phase shifts, it is more convenient to rewrite Eq. (Bl) as 


tan(5o = —Ak — Bk^ + 0{k^), 


(B2) 


where B = \A^rQ. More generally, kcotdo in Eq. (Bl) can be expanded in even powers of 


k, and tan(5o in Eq. (B2) in odd powers of k 
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In the presence of the long-range interaction —Cn/r'^ the hrst “anomalous” term in the 
expansion (B2) is proportional to [20] • This term can be calculated in the Born ap¬ 
proximation according to the prescription given by Landau and Lifshitz [20] • Consider the 


case n = 6 (van der Waals interaction). The corresponding correction A/ to the scattering 
amplitude is 


A/(g) = 2mC^q^ 


sin^ 




JqRc 

where q = 2/csin6*/2 is the momentum transfer, 6 is the scattering angle, and Rc is a cut¬ 
off radius similar to that introduced in Eq. ([^. Integrating several times by parts and 
expanding the result in powers of kRc at small kR^ gives 


A/(g) 2mC'6 




Tiq^ 


Expanding this amplitude in partial waves, we obtain for the correction to the s-wave phase 
shift; 


A5g — Cik C2k^ T 


where 


Cl = 


2mCQ 
3i?3 ’ 


C 2 = 


2TrmCek^ 

15 

2mCQ 

Rr 


(B3) 


The hrst two terms in the Eq. (B3) expression are of the same type as those in the effective- 


range expansion (B2), while the last term is “anomalous”, caused by the power-law behavior 


of the potential. Therefore, the modihed effective-range expansion can be written as 

tan(5o = —A k — B k^ -\ -h O(k^), 

15 


(B4) 


where we have introduced the new parameters A', B' > 0 to emphasize that they are different 


from those in Eq. (B2). 

One might ask if this derivation is rigorous enough because of the use of the Born ap¬ 


proximation. In fact, the expansion in Eq. (|B4|) can be derived from a more rigorous mod¬ 
ihed ehective-range theory 


for the —Cn/r” potential. This theory shows that the hrst 
“anomalous” correction to the ehective-range expansion of the phase shift can be obtained 
from the analytical continuation of the integral [2B] . 


A6l = nmCnk 


n—2 


X 


n—1 


-dx, 


(B5) 
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which converges for L > (n — 3)/2, to any physical value of L, e.g., L = 0 for the s-wave 
scattering. [In (B5), is the Bessel function.] For n = 6, one obtains (see also Refs. PUIITT] ) 


A^r, = 


Q'lrmCk'^ 


(B6) 


'■ (2L + 6)(2L + 3)(2L+ 1)(2L- l)(2L-3)’ 

which means that the lowest order correction in the effective range expansion for L = 0 is 


A5n = 


2TimCQk* 

15 ^ 


in agreement with Eq. (B3). 

For a short-range potential, the low-energy behavior of the higher partial wave phase 
shifts is 6l oc This means that for 2L-I-1 > 4, i.e., L >2, the “anomalous” correction 


(B6) is, in fact, the leading term in the low-fc expansion. This explains the behavior of the 
d-wave phase shifts seen in Figs. and when the van der Waals interaction is included. 
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